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Riassunto In questo seminario presentiamo un risultato ottenuto in collabo- 
razione con M.Manfredini e A.Sarti. Stabiliamo una relazione formale fra un mod- 
ello di oscillatori neurali nella corteccia visiva, ed il funzionale di Munford e Shah. 
Infatti le connessioni fra le colonne corticali inducono una metrica riemanniana, e la 
fase degli oscillatori evolve in questo spazio, attraverso un’equazione alle differenze, 
che possiamo pensare definita su una griglia di dimensione e. Il funzionale di Eu- 
lero Lagrange associato all’equazione T-converge, quando la dimensione della griglia 
tende a 0, al funzionale di Mumford e Shah nella metrica riemanniana. Analoga- 
mente, le soluzioni dell’ equazione discreta di fase convergono ad una soluzione BV 
che puo’ essere considerata il flusso associato al funzionale di Mumford e Shah. 
In questo modo viene quindi fornita una motivazione biologica a questo famoso 
funzionale. 


Abstract In this seminary I present a joint work with M.Manfredini e A. Sarti. 
Aim of this study is to provide a formal link between a oscillatory neural model, 
whose phase is represented by a difference equation, and the Mumford and Shah 
functional. A Riemannian metric is induced by the pattern of neural connections, 
and in this setting the difference equation is studied. Its Euler Lagrange operator 
T-converges as the dimension of the grid tends to 0, to the Mumford-Shah func- 
tional in the same Riemannian space. Correspondingly the solutions of the phase 
equation converge to a BV function, which is interpreted as the flow associated 
to the Mumford and Shah. In this way we provide a biological motivation to this 
celebrated functional. 
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1. INTRODUZIONE x 


Questo seminario è la presentazione di alcuni risultati ottenuti in collaborazione 
con Sarti e Manfredini, e contenuti in [34], dove vengono presentati vari modelli di 
visione, e posto il problema della relazione fra questi, e [13], dove viene data una 
prima risposta analitica ad alcune delle domande ivi contenute. Abbiamo studiato 
in questi lavori la relazione fra diversi modelli di ”binding” o ”perceptual group- 
ing”, che studiano processi di riconoscemento delle immagini da parte del nostro 
sistema visivo. Infatti lungo la traiettoria che va dall’oggetto fisico all’osservatore, 
le radiazioni luminose sono completamente indipendenti le une dalle altre, la retina 
è a sua volta un mosaico di cellule indipendenti. Eppure il nostro sistema visivo è 
in grado di elaborare gli impulsi che riceve, e ricostruire l’oggetto che vede. Inoltre, 
in presenza di oggetti sovrapposti è anche in grado di ricostruire bordi nascosti (il 
fenomeno e’ noto come completamento amodale). Per spiegare questi fenomeni sono 
stati introdotti sia modelli microscopici per la determinazione delle funzionalità 
biologiche [39], sia modelli macroscopici dedotti da esperimenti di fenomenologia 
percettiva (si veda [22]). 


1.1. Modelli neurali. Shuster and Wagner [38, 39] hanno sopposto che i neuroni 
nella corteccia siano raggruppati in strutture complesse, dette colonne corticali, 
ciascuna delle quali può essere interpretata come oscillatore. La fase degli oscillatori 
evolve attraverso un’equazione alle diffenrenze, e le zone semanticamente omogenee 
nell'immagine vengono codificate attraverso oscillatori che oscillano in fase (si veda 
[21]). Poiché la funzione u rappresenta la fase, ha valori in R/27 e la funzione @ è 
periodica del medesimo periodo. Inoltre il rapporto incrementale è definito a meno 
di 27 i ° 


(1) Dtu(z) := PET) EE), 


Shuster and Wagner hanno altresi proposto di convolvere con una Gaussiana che 
esprime la probabilità che un oscillatore si connesso con un altro. Ottengono quindi 
l'equazione seguente: 


6% p-e(3 
(a) dute) = f exp(DÉSED7* (Gi (Del) ) 

per una opportuna funzione @x|(z). La matrice g;;, è semidefinita positiva e de- 
scrive le connessioni fra le colonne corticali [19]. In questo primo lavoro abbiamo 
supposto che la metrica fosse riemanniana secondo il modello di [35], altri autori 
[30, 31, 32] hanno invece proposto modelli di metriche subriemanniane. 


1.2. Modelli di tipo curvatura o variazionali. E’ stato parallelamente svilup- 
pato un vasto programma di ricerca per introdurre e studiare la segmentazione di 
immagini. Assegnata un’immagine I, cioè una funzione definita su un aperto Q si 
tratta di riconoscere sottinsiemi {; che definiscono i diversi oggetti. I modelli più 
diffusi in questo senso sono il modello di Mumford and Shah [28], e i modelli di 
moto per curvatura, in uno spazio euclideo o Riemanniano. 

1.2.1 La segmentazione di immagini attraverso evoluzione per curvatura è stata 
diffusamente studiata in spazi euclidei e riemanniani. Una curva in R? puo’ essere 
identificata con un insieme di livello di una funzione u di due variabili, e a valori 
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reali. L'equazione che esprime il moto di ciascuno dei suoi livelli per curvatura in 
una varietà riemanniana si scrive 
i z o-4u) 

(3) 9;u = |Dulsd; ( ID, 
dove g'' è l’inversa della metrica. Questo approccio, è stato inrodotto da Osher 
and Sethian [29], e sviluppato da Malladi, Sethian and Vemuri [27], Caselles, Catte, 
Coll, and Dibos in [9] per segmentare immagini in spazi euclideo e e R. Kimmel, R. 
Malladi and N. Sochen in [23, 40, 41] per segmentazione in spazi riemanniani. In 
questi lavori il riconoscimento delle immagini viene effettuato evolvendo una singola 
curva nel piano bidimensionale dell’immagine. Recentemente è stato presentato da 
[35] un nuovo modello di evoluzione in cui tutti i livelli di u vengono fatti evolvere 
allo stesso tempo. La soluzione tende verso una funzione u costante a tratti con 
salti geodesici e ciascuna regione in cui u è costante individua un singolo oggetto, 
mentre i salti costituiscono i bordi degli oggetti. L’evoluzione tridimensinale di 
tutta la superficie risulta pi efficiente rispetto ai modelli precedenti, poiché, in 
presenza di immagini complesse, permette di riconoscere quali oggetti si trovano 
sopra, e quali sotto (fenomeno della saliency), e anche di ricostruire bordi nascosti 
di oggetti sottostanti (completamento modale e amodale). 

1.2.2 Il funzionale forse pit noto per il riconoscimento di immagini è il funzionale 
di Mumford e Shah, introdotto per la segmentazione di immagini in [28] 


(4) MS(u) = J |Vul? + dH"-1(S,) + J jr gli 


Di questi tre termini il primo richiede che la soluzione sia regolare fuori dall’insieme 
dei salti, il secondo impone la minimizzazione della lunghezza dei salti stessi, men- 
tre il terzo termine impone che la soluzione sia un’ approssimazione in norma L? 
dell’ immagine assegnata uo. L'operatore è definito rispetto alla metrica euclidea, 
ma come vedremo è possibile dare una formulazione analoga dell'operatore, anche 
rispetto ad una metrica riemanniana. L’esistenza di un minimo del funzionale è 
stata provata da De Giorgi, Carriero, Leaci in [18], sono studiate le sue proprietà 
di inferiore semicontinuità Anbrosio in [1] nel caso isotropo, Fonseca Fusco [20] e 
Trombetti [42] nel caso riemanniano. 


1.3. Relazione fra questi diversi modelli. Questi due modelli sono legati da 
alcune relazioni naturali. Nello studio dei minimi locali del funzionale di Mumford 
and Shah, e.le proprietà geometriche del suo insieme di salti, possiamo supporre 
che il terzo termine non sia presente. In questo caso si prova (si veda [3]) che 
i minimi locali sono funzioni costanti a tratti. Il loro insieme di discontinuità è 
stato completamente classificato soltanto in dimensione 2, ma anche in dimensione 
maggiore si sa che, se è di classe C!, allora ha curvatura nulla. Le principali 
proprietà dell'insieme del salti sono state provate da Ambrosio, Fusco, Pallara [3], 
Bonnet [8], David, in [16, 15, 17]. Si tratta quindi della medesima soluzione che si 
determina con l’evoluzione per curvatura, sia nel caso euclideo, sia riemanniano. 
Nel lavoro con Sarti e Manfredini [34], che presento qui abbiamo provato una 
prima relazione fra i modelli neuronali e quelli variazionali. La prova è organizzata 
come segue: per ogni e > 0 indichiamo È, il funzionale di Eulero Lagrange relativo 
all’equazione (2). Usando la tecnica dello slicing proviamo che questa famiglia 
di funzionali T- converge all'operatore di Mumford e Shah, con una opportuna 
metrica riemanniana. Inoltre fissato e > 0, ed un dato iniziale L? e costante a tratti, 
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il problema di Cauchy relativo all’equazione (2) ha una soluzione wu. € C([0, c0[, L?), 
e con valori costanti a tratti. Fissato t > 0, la famiglia u.(t) converge ad una 
funzione u. Rimane ancora aperto il problema di qual’è l’equazione verificata dai 
salti di u durante l’evoluzione, e se u può essere interpretata come flusso relativo al 
Munford e Shah. Cosa che invece è stata dimostrata in dimensione 1, e in assenza 
di metrica. 

La tecnica che utilizziamo è stata introdotta da Gobbino [24, 25, 26] per provare 
una congettura forumalata da De Giorgi. Egli ha provato la convergenza di una 
famiglia di funzionali alle differenze finite T convergeva al funzionale di Mumford 
and Shah, con metrica costante, e sotto ipotesi diverse dalle nostre. Chambolle [11] 
e Chambolle, DalMaso [12] hanno dato un’implementazione numerica di questa 
tecnica. La stessa tecnica è stata diffusamente estesa da Braides et alii in [4], [5, 6], 
[7] per studiare la relazione fra espressione discreta dell'energia dei mezzi elastici, 
e la sua controparte continua. 


2. DAL MODELLO DISCRETO AL MODELLO CONTINUO 
2.1. Equazione alle differenze. Studiamo dapprima l’esistenza di soluzioni per 
l'equazione (2). _ 
La matrice g;; che vi compare è continua in R”, e definisce una metrica rimanni- 


ana. La funzione @.je|, è una funzione continua, ottenuta da un opportuno rescaling 
della funzione seno: 
TT lsin(ez) if|z| < 2e7!/? 
5 :|--,-|]+R ERGE “ ; 
(5) di € ni d(2) se if |z| > 2e71/2 + ef 
Inoltre $.(z) è estesa su tutto R per periodicità. 
Poiché l’equazione è discreta, le sue soluzioni si cercano nello spazio 
PC? = {ue L?(R"): u constante sul cubo z + [0,e]" Vz e Z"}. 
Indicata 4. una primitiva di @.(2), 
patina if le] < 26/2, 
si sie 1 (sin?(4) +1+ oe) if |21>2e1/2+e4° 
il funzionale 


(6) F= f, MET f,, 3504 (DEu(2)) de), 
è definito in PC?. Inoltre 


Proposition 2.1. Sia e > 0 fissato. 
e Per ogni u € PC? il gradiente di È. in u è dato da: 


(VE) (2) = — f espl lE 07 (05 (0 (Dea) ) a 


e VÈ. è una funzione lipschtziana in PC?. 


Il problema di Cauchy associato all’equazione (2) risulta quindi: 


u(t)=-VÉF.(u.(t)) t>0, 
© { 


(0) = oe 
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dove uo € PC?. La teoria standard delle ODE fornisce una soluzione u. € 
C*([0, +00[, PC?), che dipende con continuità dal dato iniziale. 


2.2. Funzioni a variazione limitata e I convergenza. Studiamo ora la con- 
vergenza della famiglia F. di funzionali che abbiamo definito, e della famiglia 
di soluzioni u. dell'equazione differenziale. Richiamiamo preliminarmente alcune 
definizioni e proprietà di queste funzioni. (Si veda anche [14]) 

Sia 2 CR” un aperto. Indichiamo BV(0) l’ insieme delle funzioni a variazione 
limitata. E’ noto che la derivata distribuzionale Du di una funzione in questo spazio 
è una misura di Radon, che ammetta la rappresentazione seguente: 

Du = D°u + Diu + D°u, 
dove D°u è assolutamente continuo rispetto alla misura di Lebesgue e 
Diu= (u*(y) — u(y))vuH"[S(u) 
è la parte di salto, mentre D°u è la parte cantoriana della misura Du. 
Definition 2.1. Sia N C R un aperto, e sia u € BV(Q). diciamo che u è una 
funzione speciale a variazione limitata, e indichiamo SBV(Q) se D°u = 0. Diciamo 
che u € SBVioc(0) se u € SBVioc(A) per ogni A CCL (si veda [2]). 


Ricordo per completezza anche la definizione di I-convergenza di funzionali 


Definition 2.2. Se (X,d) è uno spazio metrico, una famiglia F; : X + R di 
funzionali I converge a F per j + 00 se le due condizioni seguenti sono soddisfatte: 
e Per ogni u in X e per ogni successione (u;) convergente a u in X, 


F(u) < liminf F;(u;) 
i 
e per ogni u € X esiste una successione (u;) convergente a u in X tale che 
F(u) > lim sup F;(u;). 
j 
Per la particolare espressione dei funzionali f., avremo bisogno di una opportuna 


modifica del classico teorema di compattezza negli spazi BV. Infatti lavoreremo 
con un a famiglia (w.) di funzioni tali che la quantità 


@® N(u) =  luldz+ fe! |. 1Dgu.(a)idede 


sia limitata 

Theorem 2.1. Sia (u.) una famiglia di funzioni in Li,, tali che N(u.) è limitata. 
Allora esiste una successione e; convergente a 0 e una funzione u in BVioc tale che 
Ue; converge a u în Li. 


La prova si può fare verificando che la famiglia 
ui(a) = f n(eru.(e + «Md 
è limitata in BV(Q), e quindi relativamente compatta. D'altra parte 
fu = ua) < ce 


e quindi anche la famiglia assegnata ammette una sottosuccessione convergente. 
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2.3. convergenza di funzionali di dimensione 1. La prova della convergenza 
della famiglia di funzionali F; che presentiamo è basata sul metodo delle sezioni. Si 
prova preliminarmente un risultato di conergenza in dimensione 1. Poi, sfruttando 
il fatto che i funzionali sono in forma integrale, si estende il risultato in dimensione 
alta, applicando un’opportuna variante del teorema di riduzione, che vale per le 
funzioni in esame. La tecnica che usiamo combina quelle di [24] e [4]. 

In dimensione 1 l’operatore diviene 


(O I Fogg(1) = f 9(2)04g (Diu(2)) de 


dove I C R, g è una funzione continua e strettamente positiva. In dimensione 1, 
fissata £ positiva si ha 


J 9() x] (Diu(2)) de = J 9(201g) (Difu(x)) de, 


e analogamente si può ragionare per £ < 0. Quindi supporremo sempre che € = 1 
in questo contesto. 
Ricordiamo che 


Theorem 2.2. Se 4 : [0,+00[+ R è una furizione convessa non decrescente 
lim:s4o ell = +00 e % : [0,+00[+ R è una funzione continua, subadditiva e 
tale che inf $ > 0, allora il funzionale MSp,w,g(u, I), definito 


(10) 
Sro(a)e(IVu(2))ez + Sinsqy A(2V(u* — dr w € SBV() 
MSo,4,9(U, I) = 


+00 si altrimenti 


è inferiormente semicontinuo in Li,.- 


i Theorem 2.3. Indichiamo 
2° sin?(2) 
g()= 7: v@)==774+1 
Allora per ogni successione u; + u în Li, si ha 
lim ai Fe;,1,g(0;,R) > MSpw,g(u, R). 


lim sup F,,ig(u,R)< MSo,y,g(u, R). 


j++t00 
Plim Fego(u;,R) = MSp,y,o(u,R). 
diamo un accenno di prova della prima affermazione. Possiamo supporre, ragio- 


nando come [4] p. 82, ed eventualmente modificando la successione u., che per ogni 
n > 0 sia verificata la relazione 


Fi,(u,1) +n2 6; D_ g(ke;)0x; (D'u(ke;)), 
keJ; 
dove abbiamo indicato 
Jj = {k eZ: (ejk, e;(kK + 1)) E I}. 
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Poiché J;j è finito, possiamo indicare 
J; = {K,---Ky,} 
e denotiamo 
u;((k +1)6;) — u;(ke;)) mod2r =) 
I = {Ke J;: Lt +16) — u;(ke;)) mod2r] s(( i) a s( ci < 26; sd" Jj=Jj\J}. 
Poi definiremo v; come segue: 


(£ -_ k)u;(e;(k + 1)) + (( + 1) — 1 )u;(k6;) in €;(k,k la 1), k E Jj 


u;(ke;) în e;(k,k + 1), ke Jî 
u;(k66;) if t< koe;. 
u;((kf, + 1)6;) if t> (ky, +1);. 


Con queste notazioni la stima di F. ; diviene 


Fi;(u,1)+n2 6 YO g(ke;)px, (DI u(ke;)) > 
keJj 


2 DO Gig(ke;)pa (DSu(k6;)) + DD iglte)p,, (LETT R)) , 
keJ} keJ} 3 


Le ejg(ke;)p(D'u(ke;))+ YO 9(Ke;)Y (4; ((k+1)€;)-u;(ke;)) 


= —1/2 8 
“6; #(8) keJ} keJ? 


= inf Le; (8) 9(t).0(|v;|)dt + > 9U(lv} — v7 |). - 


ss! (8) Jr S(v;)0I 
Per il teorema precedente M S è inferiormente semicontinuo, e quindi si ottiene 
lim inf F,(u;,I)> MS(u,I)-n. 
j++00 


L’ arbitrarietà di 7 > 0 fornisce la prima affermazione nel caso di / intervallo 
limitato. L'affermazione rimane vera anche su tutto R come si prova approssimando 
questo insieme dall’interno. 


2.4. Il metodo delle sezioni: convergenza di F. in dimensione n. Il caso di 
dimensione n segue da quello di dimensione uno, con il metodo di slicing, diffusa- 
mente applicato a problemi di convergenza per funzioni BV. Si veda per esempio 
Gobbino [24], and Braides [4]. Riportiamo qualche dettaglio della prova, per poter 
definire l'operatore di MS rispetto ad una metrica riemanniana. 7 

Poiché abbiamo necessità di calcolare il rapporto incrementale nella direzione 
£, rappresentiamo lo spazio come prodotto cartesiano di una retta parallela a £, e 
dell’ortogonale. Ora infatti £ € R"\{0} e sia < £ >1= {z € R° :< €£,2 >= 0} lo 
spazio ortogonale a £. Per ogni y €< € >+ consideriamo la funzione Ugy definita 

Ugy = u(y +66), te R. 


Con queste notazioni l'operatore f. definito in (6) diviene 


(11) Fi R") = [EE [Fans ey, R)dudk, 


Il 
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dove l’operatore F. €,9;; (Uey, R) è stato definito nella sezione precedente. In questo 
modo il funzionale Fi 6,9; si rappresenta in termini delle sue sezioni unidimensionali. 
Analogamente anche il funzionale di MS si rappresenta in termini delle sue sezioni. 
Introduciamo ora la metrica 
—1elE;E; 
e @ 
(12) YijPiPj = [rif p, € >)? dé 


e definiamo l’operatore di Mumford e Shah in dimensione n: 
(13) = MSsyg(u,R")= f sd; - / W(lut — u7|)Ny(x, vu)dH"-!, 
S(u) 


dove : 


(14) PASTA / . e MES 3,0] <v,€ > |de. 


Theorem 2.4. Per ogni funzinoe u € Lì .(R") si ha 


loc 
MSpyg(R")= f eMÈS(/ MS6g4 (Ue R)dy)d, 


n 1? 
dove MSj,y,9 è definita in (13). 


<E>+ 


Le proprietà del funzionale È. si deducono ora dai risultati unidimensionali e 
dalla formula di rappresentazione contenuta nel Teorema 2.4. 


Theorem 2.5. La famiglia di funzionali (É.) definita in (11) T- converge a 
MSo,6,9 în bi 

2.5. Convergenza delle soluzioni dell'equazione di evoluzione. Proviamo in 
questa sezione che, fissata una funzione uo, questa puèssere approssimata da una 
successione u. di funzioni tali che u. € PC? per ogni e. Per ogni e il problema di 
Cauchy introdotto nel paragrafo 2.1 individua una funzione (ue(t)). Qui proviamo 
che (u.(#)) ha una sottosuccessione convergente ad una funzione u(t) tale che u(0) = 
uo. I limiti possibili di successioni di questo tipo sono i candidati ad essere il flusso 
gradiente del funzinale di Mumford e Shah con dato iniziale uo. 

Definiamo lo spazio 


X = {u € SBVio:(R",R/[0,27)) : MS(u) < +00). 


Poiché il funzionale È. è definito per funzioni costanti a tratti, introduciamo una 
discretizzazione della funzione uo. 


Proposition 2.2. (si veda [25]). Esiste una famiglia (uoe) € PC?(R",R/(0,27]) 
tale che 
uo duo in L}.(R"); 


lim F,(uoe) = MS(uo) 
€+0 


sup{F} (uoe)} < +00. 
€>0 


Studiamo ora il comportamento della famiglia (w.). 
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Lemma 2.1. La famiglia (u.) determinata nella sezione 2.1 soddisfa 


(15) sup (/ fuolde + f et |Dju,(2)|dzd£) < +00 
e€[0,1] \J9 R” lo 


ed esiste una successione (ex) convergente a 0 tale che (ue,) è relativamente com- 
patta in C((0,+00[; L}..(R")) e ha pertanto un limite u € C([0,+00[; Lî.c(R")) tale 
che u(t) € BV(Q) per ogni t. 


Notiamo preliminarmente che il funzionale È. è decrescente lungo il moto ue. 
Infatti 


(16) Sf (ue(0)) = (VÉ (0: (0)) UL) 2) = 


= —|lu(t)]]z2 e») = — VE (ve(t))]l22 8") ! 
Vista l’ espressione della funzione g., è naturale definire Di*u.(7) = Dju.(x) se 
|Déue(2)| > (el) e Di'*u.(x) = 0. Inoltre chiamiamo 


(17) I = {2 e R*IDf+u.(2) #0). 


Per ogni N CC R", poiché u. assume valori in [-, 7], si ha 


I ={ eis Li |Déue(x)|dedf < 2e f et ls 1A qqidi < 
n toa R” rtna €lE 


< 2e / em ls ele (Djue)drdf < cF.(u.). 
n TÈN9 
Inoltre poiché g;; è uniformemente ellittica, 


F.(w.) > :[, Lo e lp.|e(Djue(2))drdt > cf. ta ellp(Dju.(2)) dedé 


R"\I% 


per definizione di £. ‘Ne viene che 


ri ell J, |Djus(x)|drdf = I + / e J |Déu(x)|drdt < 
È Q pr QI 


< cF;(u) + (A el6! (Lu 


«€ 


< cF(ue) + c|Q, 


(Diu, (ode + 1) de < 


per una costante opportuna c. 

Pertanto (15) è dimostrata. Da teorema 2.1, segue che per ogni # > 0, la famiglia 
(u.(t)) è relativamente compatta in L?(R"). 

Inoltre, 


ta t2 } 
1 
Ilue(t1)-ue(t2)]|22(R") </ luce(t)]|zz0en)dt < (/ lu (ica) ltit2]? < 
ti 


ti 


1 
n 2 1 1 
< (Fe(o(8)))” li — tal? <clti — tal 


per ogni e > 0, e la conclusione segue dal teorema di Ascoli-Arzelà. 


Le principali proprietà della funzione limite sono le seguenti: 
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Theorem 2.6. Sia ue soluzione del problema (7), e u il suo limite, allora: u € 
C([0, +00; Lioe(R")) u(0) = uo; MS(u(t)) £ MS(uo); la funzione u = u(7,t) è 
soluzione distribuzionale ]0, +00[xR" dell ‘equazione 


du 
n” D(g:; Vu) 


dove D è la derivata distribuzionale fatta rispetto ad x, Vu è il gradiente approssi- 
mato di u. 
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